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Abstract 

For a closed system to equilibrate from a given initial condition there must exist an 
equilibrium state with the energy equal to the initial one. Equilibrium states of a 
strongly coupled gauge theory with a gravitational holographic dual are represented 
by black holes. We study the spectrum of black holes in Pilch-Warner geometry. These 
black holes are holographically dual to equilibrium states of strongly coupled SU(N) 
Af = 2* gauge theory plasma on S 3 in the planar limit. We find that there is no energy 
gap in the black hole spectrum. Thus, there is a priory no obstruction for equilibration 
of arbitrary low-energy states in the theory via a small black hole gravitational collapse. 
The latter is contrasted with phenomenological examples of holography with dual four¬ 
dimensional CFTs having non-equal central charges in the stress-energy tensor trace 
anomaly. 
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1 Introduction and summary 

Consider an interacting system in a finite volume. Suppose that the theory is gapless — 
there are arbitrarily low-energy excitations. If a generic state in a theory equilibrates, 
there can not be a gap in the spectrum of equilibrium states in the theory. This 
obvious statement has a profound implication for strongly coupled gauge theories with 
an asymptotically AdS gravitational dual [1], In a holographic dual the equilibrium 
states are realized by black holes [2], Thus, if it is possible to prepare an arbitrary 
low-energy initial configurations in a holographic dual with a gapped spectrum of black 
holes, such states of the boundary gauge theory will never equilibrate. Correspondingly, 
the asymptotically AdS dual is guaranteed to be stable against gravitational collapse 
for sufficiently small amplitude of the perturbations. Examples of this type would 
violate ergodicity from the held theory perspective. 
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In this paper we show that while it is possible to realize above scenario in a phe¬ 
nomenological (bottom-up) holographic example — the Einstein-Gauss-Bonnet (EGB) 
gravity with a negative cosmological constant, it does not occur in a specific model 
of gauge theory/supergravity correspondence we consider — the holographic dual¬ 
ity between Af = 2* SU(N) gauge theory and the gravitational Pilch-Warner (PW) 
flow [3-5]. 

From the gauge theory perspective, SU (. N) Af = 2* gauge theory is obtained from 
the parent Af = 4 SYM by giving a mass to Af = 2 hypermultiplet in the adjoint 
representation. In R 3 ' 1 space-time, the low-energy effective action of the theory can 
be computed exactly [6]. The theory has quantum Coulomb branch vacua A4c, pa¬ 
rameterized by the expectation values of the complex scalar in the Af = 2 vector 
multiplet, taking values in the Cartan subalgebra of the gauge group, 

$ = diag(ai,a 2 , • • • ,a N ), ^^ = 0, (1.1) 

i 

resulting in complex dimension of the moduli space 


dim c A4c = N — 1. 


( 1 . 2 ) 


In the large -N limit, and for strong ’t Hooft coupling, the holographic duality reduces 
to the correspondence between the gauge theory and type lib supergravity. Since 
supergravities have finite number of light modes, one should not expect to see the full 
moduli space of vacua in Af = 2 examples of gauge/gravity correspondence. This is 
indeed what is happening: the PW flow localizes on a semi-circle distribution of (1.1) 
with a linear number density [4], 


Im(aj) = 0, 
8n 


&i £ [ — Q-o, Clo] , 


a 0 — 


m 9 ym n 
47T 2 


p(a) = 


m2 9YM 


ra 0 


da p(a) = N , 


(1.3) 


r —CLQ 


where m is the hypermultiplet mass. This holographic localization can be deduced 
entirely from the field theory perspective [7], using the ^-supersymmetric localization 
techniques [8]. To summarize, Af = 2* holography is a well-understood nontrivial 
example of gauge/gravity correspondence that passes a number of highly nontrivial 
tests [4,7,9]. 

We would like to compactify the background space of the Af = 2* strongly coupled 
gauge theory on S 3 of radius f — in a dual gravitational picture we prescribe the bound¬ 
ary condition for the non-normalizable component of the metric in PW effective action 
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to be that of R x S 3 . This is in addition to specifying non-normalizable components 
(corresponding to m in (1.3)) for the two PW scalars, dual to the mass deformation op¬ 
erators of dimensions A = 2 and A = 3 of the gauge theory hypermultiplet mass term. 
Thus, we produced a holographic example of a strongly interacting system in a finite 
volume. The single dimensionless parameter 1 , so far, is ml. We proceed to construct 
regular solutions of the PW effective gravitational action with the prescribes bound¬ 
ary condition, interpreting them as vacua of Ahcompactihed strongly coupled J\f = 2* 
gauge theory. Using the standard holographic renormalization technique 2 we compute 
the vacuum energy of the theory as a function of m£, E vacuum = E vacuum (m£). We do 
not verify in this work whether described 5' 3 -compactihcations preserve any supersym¬ 
metry; thus, it is important to check the stability of the vacuum solutions. Previously, 
careful analysis of the ^-compactihed PW holographic flows of [11] pointed to the dis¬ 
crepancy in the free energy of the solutions, compared with the localization prediction 
in [7]. This discrepancy was resolved by identifying a larger truncation [9] (BEFP) 3 , 
where it was pointed out that preservation of the ^-supersymmetry necessitates turn¬ 
ing on additional bulk scalar fields. Stability of the PW embedding inside BEFP was 
discussed in [12], We verify here that S' 3 -compactihed PW vacua are stable within 
BEFP truncation. Having constructed vacuum solutions, we move to the discussion 
of the black hole spectrum. We construct regular Schwarzschild black hole solutions 
in PW effective action, and compute 5E = SE(m£, £bh / L) = E — E vacuum (ml). We 
argue that there is no obstruction of initializing arbitrary low-energy excitations over 
the vacuum. Thus, one would expect no gap in the energy spectrum of PW black hole 
solutions, realizing equilibrium configurations of the strongly coupled Af = 2* gauge 
theory in the planar limit. Indeed, we find strong numerical evidence that 


lim 5E(m£, £ B h/L) 
£bh/l-^§ E vacuum (m£ 0) 


0 . 


(1.4) 


The rest of the paper is organized as follows. In the next section we discuss the 
spectrum of black holes in five-dimensional EGB gravity with a negative cosmological 
constant. These gravitational backgrounds can be interpreted as holographic duals to 
equilibrium states of strongly coupled conformal gauge theories with non-equal cen¬ 
tral charges in the stress-energy tensor trace anomaly. We show that there is a gap 

1 Af = 2* theory in Minkowski space-time has a scale associated with the Coulomb branch moduli 
distribution (1.3). Once the theory is compactified on the S 3 the moduli space is lifted. 

2 For the model in hand this was developed in [10]. 

3 Of course, BEFP can itself be consistently truncated to PW. 
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in the spectrum of black holes. However, as one imposes constraints on EGB gravity 
coming from interpreting it as an effective description of gauge theory/string theory 
correspondence, the claim about the gap becomes unreliable — higher derivative cor¬ 
rections, which are not under control, make order-one corrections to the gap. We follow 
up with the discussion in the Af = 2 * holographic example. In the section 3 we review 
the PW effective action and its embedding within a larger BEFP truncation. In sec¬ 
tion 4 we construct gravitational dual to vacuum states of Af = 2 * gauge theory on 
S 3 . Stability of the latter states within BEFP truncation is discussed in section 5. In 
section 6 we study the spectrum of black holes in PW effective action. 


2 Black hole spectrum in Einstein-Gauss-Bonnet gravity 


Effective action of a five-dimensional Einstein-Gauss-Bonnet gravity with a negative 
cosmological constant takes form: 

S= ^lL d5z ^{^ + R+ ^ L2 ( R2 - 4R ^ U + R ^ R,XUpa ))- ( 21 ) 

When interpreted in a framework of gauge theory/gravity correspondence 4 , EGB action 
( 2 . 1 ) represents a holographic dual to a putative strongly coupled conformal theory with 
non equal central charges, c ^ a, of the boundary stress-energy tensor, 

Ea = r^p X r^ px - 4+ r 2 , (2.2) 

h = V^-2v' fl '+^ 2 i 

where E 4 and I 4 are the four-dimensional Euler density and the square of the Weyl 
curvature of the CFT background space-time. The precise identification of the central 
charges is as follows: 


7 r 2 Z 3 


c = 


£ 3 


1-2 


A 


GB 


L = ^L, 


P 2 


7 t 2 L 3 


a = 


i 3 


1 - 6 - 


*GB 




4A, 


GB 


(2.3) 


The gravitational dual to the vacuum state of a CFT on a three-sphere S 3 is a global 


AdS 5 , 


ds 2 = 


L 2 (3 2 

COS 2 X 


(— dt 2 + dx 2 + sin 2 a;dr2 2 ) , xe[0, 7t/2], 


(2.4) 


^See [13] for a recent review. 
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where dVL\ is the metric of S 3 . Notice that A G b is restricted to be 


A 



(2.5) 


otherwise, there is simply no asymptotic AdS solution. Following holographic renor¬ 
malization of EGB gravity developed in [13,14], we find that the vacuum energy (the 
mass) of (2.4), or the Casimir energy from the boundary CFT perspective, is 


E 


vacuum 


3 a 
4L ' 


( 2 . 6 ) 


Black holes (equilibrium configurations of EGB CFT) are found as a regular horizon 
solutions within the metric ansatz, 


L 2 f3 2 

COS 2 X 


ds 2 = 


d X 2 

-A(x)dt 2 + „, , + sin 2 xdQ, 


A(x) 


(2.7) 


The most general solution of equations of motion obtained from (2.1) determine A(x) 
is terms of a single parameter M > 0, 


A 


1 - —— ((2X G b - /3 2 ) sin 2 x + (a\ gb (/3 2 - 2A GB )M cos 4 a; 

2\GB \ \ 

+ (2 X gb — (3 2 ) 2 cos 4 x — /d 4 (l — 4A gb ) cos(2x) 



( 2 . 8 ) 


Furthermore, using the machinery of the holographic renormalization, the energy of 
the boundary CFT is 


E 


3c f f3 2 — Q\ G b 
U^\^-2X gb 


+ 4 M 


3c 

41 


+ 4 M 


(2.9) 


It is remarkable that the regular Schwarzschild horizon in the geometry (2.7), (2.8) 
exists only provided [15,16] 


M > I 2 ^- 4 ’ 

1 {( 3 2 ~ 1 ) ( 2 / 3 2 


if A gb > 0 , 
1), if A gb < 0 . 


( 2 . 10 ) 


For positive A G b, the bound comes requiring that S 3 remains finite at the location 
of the horizon (otherwise the curvature at the horizon diverges). For negative A G b, 
violating the bound would render geometry complex (expression inside the square root 
in (2.8) would turn negative for some x G (0, 7t/2)). 
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Constraints (2.10) imply the gap in 8 E = E — E vacuum in the spectrum of EGB 
black holes, 


5E 


I E 


— e gap 


vacuum 


4(1 - /3 2 ) I 1, A gb > 0, 

|6/3 2 - 5 I X 1 -(2/3= - 1 ) 2 , \ GB < 0, 


(2.H) 


With the only restriction (2.5) on A gb, £ g a P is unbounded as A gb —>■ —oo and A gb —$■ 
5/36. 

We argue now that attempts to interpret EGB holography as an effective descrip¬ 
tion of some gauge theory/string theory correspondence make the gap claim (2.11) 
unreliable. First, causality of the holographic GB hydrodynamics requires that [17] 


7 

36 


< A GB < 


9 

100 


( 2 . 12 ) 


Additionally, it was pointed out [18] that pure EGB gravity with a negative cosmo¬ 
logical constant can not arise as a low-energy limit of a gauge theory/string theory 
correspondence — the difference of central charges (c — a)/c is bounded by A J 2 p , 
where A gap is the dimension of the lightest single particle operators with spin J > 2 
in the holographically dual conformal gauge theory. Integrating out massive J > 2 
spin states generic-ally produces new higher-curvature contributions, in addition to 
the Gauss-Bonnet term. These higher curvature corrections are as important as the 
Einstein-Hilbert term and the GB term in (2.1) when the size of a black hole becomes 
of order A gbL. The latter is true even as A gb <C 1, as the Ricci scalar evaluated on 
the horizon of ~ A gbL size black hole (2.8) diverges as A—. 


3 PW/BEFP effective actions 


We begin with description of the PW effective action [3]. The action of the effective 
five-dimensional supergravity including the scalars a and y (dual to mass terms for 
the bosonic and fermionic components of the hypermultiplet respectively) is given by 


S = 


-PW 


[ d,ev=g £/ 

Jm 5 

1 d^^g[\R-3(da) 2 - (<9y) 2 -V\ , 


4vrG, 


5 J 


(3.1) 
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where the potential 5 


V = — 
16 


1 fdW 


3 V da 


+ 


dW 

dx 


\ w2 < 


(3.2) 


is a function of a and x, and is determined by the superpotential 

W = -e~ 2a - ie 4 " cosh(2x). 

In our conventions, the five-dimensional Newton’s constant is 

r _ G w _ 4tt 
5 2 5 V 0 I 55 N 2 

Supersymmetric vacuum of AT = 2 * gauge theory in Minkowski space-time is given by 
dsl = e 2A (—dt 2 + dx 2 ) + dr 2 , p = p(r) = e a ^ , % = %(r), (3.5) 


(3.3) 


(3.4) 


with 


e A = 


kp 2 


p 6 = cosh(2y) + sinh 2 (2y) In 


f = (3.6) 

cosh(x) dr 3 


sinh(2x) ’ 

where the single integration constant k is related to the hypermultiplet mass m ac¬ 
cording to [4] 

k = mL = 2m . (3.7) 

The BEFP effective action [9] is given by 


S , 


BEFP 


/ Gbefp 


An Gk 


[ dev^g 

\r- 12< S f-4 < a h 2 _vl 

)Ms 

V 2 (1 - X 2 ) 2 


(3.8) 


with the potential 


V= - 


1 10^ + ^ 8 (Ad ) 2 + (X 2 ) 2 

V 4 1-X 2 ( 1-X 2 ) 2 


(3.9) 


where X = (Ad, A" 2 , X 3 , X 4l X 5 ) are hve of the scalars and r/ is the sixth. The symmetry 
of the action reflects the symmetries of the dual gauge theory: the two scalars (Ah, Ad) 

5 We set the five-dimensional supergravity coupling to one. This corresponds to setting the radius 
L of the five-dimensional sphere in the undeformed metric to 2. 



form a doublet under the U( 1)r part of the gauge group, while (X 3 , X 4 , X 5 ) form a 
triplet under SU(2) v and 77 is neutral. The PW effective action is recovered as a 
consistent truncation of (3.8) with 

X 2 = X 3 = X 4 = X 5 = 0, (3.10) 

provided we identify the remaining BEFP scalars (rj,Xi) with the PW scalars (cqy) 
as follows 

e a = V , cosh 2x = | ^ . (3.11) 

Note that once m ^ 0 (correspondingly Xi ^ 0), the U( 1)r symmetry is explicitly 
broken; on the contrary, SU(2)y remains unbroken in truncation to PW. 

4 Holographic duals to Af = 2* vacuum states on S' 3 

We derive bulk equations of motion and specify boundary conditions representing grav¬ 
itational dual to vacuum states of strongly coupled Af = 2 * gauge theory on S 3 . We 
assume that the vacua are »S'0(4)-invariant. We argue that there is no obstruction of 
exciting these vacua by arbitrarily small perturbations of the bulk scalar fields a and y. 
We review holographic renormalization of the theory and compute the vacuum energy. 
Next, we solve static gravitational equations perturbatively in the mass deformation 
parameter ml < 1 — this would serve as an independent check for the general 0 (m£) 
numerical solutions. We conclude with the plot representing e = E vacuum {m £)/ 

3 n 2 

Evacuum = ^vacuum( tu£ = 0 ) = ’ ( 4 - 1 ) 

as a function of mi. Interestingly, while the vacuum energy of the J\f = 4 SYM is 
positive, it is negative 6 for Af = 2* gauge theory once ml > 0.87. 

4.1 Equations of motion and the boundary conditions 

We consider the general time-dependent SO (4)-invariant ansatz for the metric and the 
scalar fields: 

ds\ = — f—Ae~ 2S (dt) 2 + + sin 2 a;(dfl 3 ) 2 ^ , (4.2) 

cos^x \ A J 

6 Prior to imposing causality constraints in EGB gravity, its vacuum energy becomes negative once 
A gb > 5/36. Vacuum energy of a different nonconformal gauge theory on S 3 was also observed to be 
negative in [19]. 
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where (dfl 3 ) 2 is a metric on a unit 7 round S 3 , and {A, 5, a, x} being functions of a 
radial coordinate x and time t. Introducing 

4+ = d x a , $ x = <9 x x, n Q = — d t a , n x = — d t \ , (4.3) 

we obtain from (3.1) the following equations of motion: 

■ the evolution equations, ' = d t , 


a = Ae d U a , x = Ac d II x , 

= (Ae-‘ n«) 4> x = +-+) , 


ru = 


(tan 3 xAe 5< I> a ) — 


tan 3 x 

fly = -77— (tan 3 xAe~ 5 $A 

x tan 3 a: v x ’ 


dV 


x 3 cos 2 x da 


-5 


(4.4) 


—6 


dV 


’ x cos 2 x dx ’ 


the spatial constraint equations, 
2 + 2 sin 2 x 


A 


Sin X COS X 


— 4 tan x 


(1 — A) — 2 sin(2:r)A 


! + |p 


$ 2 + n 2 + -& 2 + -n 2 



(4.5) 


and the moment constraint equation, 
A )t + 4sin(2:r)A 2 e _<s 


$ a n Q + -$ x n x )=o. 


(4.6) 


It is straightforward to verify that the spatial derivative of (4.6) is implied by (4.4) 
and (4.5); thus is it sufficient to impose this equation at a single point. As x —* 0 + , 
the momentum constraint implies that A(0,t) is a constant 8 , and as x —> the latter 
constraint is equivalent to the conservation of the boundary stress-energy tensor (see 
4.2 for details). 

The general non-singular solution of (4.4), (4.5) at the origin takes form 


A(t, x) — 1 + 0(x 2 ) , S(t, x) = dg (t) + 0(x 2 ), 

a(t, x) = aft(t) + 0(x 2 ) , x(t, x) = Xo(t) + 0(x 2 ). 

'We set l = 1; the t dependence can be easily recovered from dimensional analysis. 
8 In fact, the non-singularity of A(t,x) in this limit automatically solves (4.6). 
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It is completely characterized by three time-dependent functions: 


{d h o,a h 0 , X h o}- (4-8) 

At the outer boundary x = § we introduce y = cos 2 x so that we have 


A—l + y -ci, 0 + y 2 [a 2 , 0 (t) + ^-ci i0 (ci )0 + 1) + 8p 2 x + 16pi,ipi i0 (i)J In y 

+ 8pi i In 2 y \ + 0{y 3 In 3 y), 

<5 =V g c i,o + — —c 2 o + 4 --Pipit) — -Ci >0 + 2 \p\ l + 4p li0 (t)pi jl 

+ Qci,o + ic 2 i0 + 4p 2 a + 8pi,o(t)pi,i^ In y + 4p 2 ;1 In 2 y^j + C>(y 3 In 3 y), 

e“ =1 + y (pi,o(t) + pi,i In y) + y 2 ^c 2 0 + pi, 0 (t) - 3pi,iCi )0 + &p\ tl 
4 3 i /4 

— 4p 1 , 0 (i)p 1 ,i + -ci,opi,o(t) + ~plfi{t) + -<9^pi, 0 (t) + ( -pi,iCi )0 + pi,i 

- 4p 2 x + 3pi i0 (t)pi,i^ In y + ^p 2 1 In 2 y^j + 0(y 3 In 3 y), 

cosh 2y =1 + y ci >0 + y 2 ^c 2 , 0 (t) + Qci >0 + ln v'j + 0(y 3 ln 2 y), 


(4.9) 


where we explicitly indicated time-dependence, i.e., 

4^ = °, 4«,=0. (4.10) 

Asymptotic expansion (4.9) is completely characterized by two constants 9 {pi,i,ci ; o} 
and three time-dependent functions 

{ fl 2,0 , Pi,o ; c 2,o} , (4-11) 


constraint by (4.6) to satisfy 

0 = ^ 2,0 - 8p 2 j0 (t) - 16pi, 0 (£)pi,i - §<*,(*)) • (4-12) 

9 Prescribing time dependence to these coefficients amounts to study quantum quenches in J\f = 2* 
gauge theory [20]. 
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The non-normalizable coefficients p u and c^o are related to the mass deformation 
parameters of the dual gauge theory. Following [21], the precise relation can be es¬ 
tablished by matching the asymptotics (4.9) with the supersymmetric PW RG flow 
(3.6), 


{pi,i> cyo} 


PW 


1 1 


48 8 


= k o 


1 1 


12 2 


(4.13) 


A specific relation between the non-normalizable coefficients of the bulk scalars e a and 
cosh2y, i.e., 


Ci,o — 6pi,i, 


(4.14) 


realizes Af = 2 supersymmetry of the boundary gauge theory in the UV. As in [21], it 
is possible to study the theory with explicitly broken supersymmetry, i.e., 

1 111 

Pu = (■ m b L) 2 ^ - X Ci,0 = - X - (m f L) 2 , (4.15) 

where nib and m/ are the masses of the bosonic and the fermionic components of the 
Af = 2 hypermultiplet of the boundary gauge theory. 

A non-equilibrium state of the gauge theory can be specified with the following 
initial/boundary conditions: 


a(0, *) = <«"**(*), x(0, x) = x inU (x ), $«(0, x) = 4>r = 


da 


init 


dx 


init d X 


init 


$ x (o,x) = $ 7 * = 


dx 

and as y = cos 2 x —>■ 0, 


a imt (y) = pi,i ylny + 0 (y ), cosh (2 x tmt (y)) = 1 + y c lfi + 0 (y 2 In y ), 
n7*(P) = 0 (y) , n ™*(y) = 0 (y 3 ' 2 ), 


A(0,x) =1 + - 


cos 4 X 


sm 2 x 


eX P 1 “g 


df sin( 2 o (nr ft (0) 2 + (^ Bit (0) : 


+ 3(n:t )) 2 + 3($ i ;t)) 2 x g(x), 


g{x) = - o / tan3 ^ ex P o 


+ 3 (n^fa)) + 3 J X 

+ (*f*(0) 2 + 3 (ic*(0) 2 + 3 K mt (0) ; 

v init (0 = r(a init (0,x init (0), 


dy sin( 2 ^(n ^(v)) +(^r t (v)Y 

AV ini \ff) + 3 


cos 2 £ 


+ ( n f“K))' 


(4.16) 


n„(o,i) = n“‘(i), n x (o, x) = n*”“(i), 


(4.17) 


(4.18) 
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,5(0,x) = sm(2{)((nf‘({)) 2 + (tf‘K )) 2 + 3 (nf‘K )) 2 + 3 (*r‘«)) 2 ) , 

(4.19) 

where we explicitly solved for AL(0,a;) and 5(0, x) using constraint equations (4.5). 
Notice that while A(0,x) and 5(0, x) are free from the singularities given arbitrary 
profiles (4.16), a large amplitude initial conditions might cause A( 0 ,x) to vanish for 
some 0 < Xq < t(, i.e., AL(0,a;o) = 0, — this corresponds to ’putting a black hole in 
the initial data’. Clearly, initial conditions arbitrarily small perturbed about static 
gravitational solutions without a horizon (see below) are well defined. In particular 
one can can consider perturbations with 

a init = a\ x init = x\ n£« - A n a Jx) , A -> 0, (4.20) 

where the superscript v stands for a static (vacuum) solution and A characterizes an 
overall amplitude of the perturbation with given initial profiles ir a and 7r x . 

The S'0(4)-invariant vacua of strongly coupled Af = 2* gauge theory correspond 
to static solutions of (4.4)-(4.6). To avoid unnecessary cluttering of the formulas, we 
omit the superscript v , use a radial coordinate y = cos 2 a:, and introduce 


A(t,y)=a(y), S(t,y)=d(y), = p(y ), cosh(2x(t, y)) = c(y). 

(4.21) 

We find then 

n = r " _ c ( c ') 2 , j (d_ _ A _ (y + ^ _ P 2 (c 2 ~ l)(p 6 c-4) 

' c 2 -l \a ) y(l-y) 4(1 - y)y 2 a ’ 


0 = p ''-^ + p ' L ±-d' - 


(y + 1 )rt (c 2 - i)p 9 


1 — p 6 c 


2/(l-p) 12(1 ~y)y 2 a 6 p 3 y 2 a(l~y) : 


Q = d ,_ m-y)(C) 2 _ 8(1 - y)y{p'f 
3(c 2 — 1) p 2 

, , 2n / 8(p') 2 2(c') 2 A 

0-a -(y-y)a (— + 3(c 2 _ j)) + 


(y — 2 )a + y (c 2 — l)p 8 — 8p 2 c 2 

y (i -y) % + 3 yp A ’ 

(4.22) 


where ' — 4~- The boundary conditions as y —» 0 are as in (4.9), once we neglect the 
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time dependence. At the origin, using z = 1 — y we have 


a — 1 + ^ 1 + 

d = dj + Ofz 2 ), 


3(pS) 4 


(M ((c 'y _ 1) + LtMl!) 2 + 0 (z 2 ), 


c = c$ + LpS) 2 ((cj) 2 - 1 ) (cS(pS) 6 -4 ) z+ 0 {z 2 ). 


(4.23) 


We consider geometries with J\f = 2 supersymmetry in the ultraviolet, so we impose 
the constraint (4.13). Having fixed m, the complete set of normalizable coefficients in 
the UV/IR is given by: 


{« 2,0 , P 1,0 5 c 2,0 


Po 


lh 1 

'0) U 0J 


(4.24) 


Note that the six integration constants (4.24) is exactly what is needed to uniquely fix 
a solution of a coupled system of two second-order and two first-order ODEs. 


4.2 Holographic renormalization and the vacuum energy 


Holographic renormalization of RG flows in PW geometry was discussed in [10]. Here 
we apply the analysis for the gravitational solutions dual to vacua of J\f = 2* gauge 
theory on S 3 . 

The gravitational action (3.1) evaluated on a static solution (4.22) diverges - 
this divergence is a gravitational reflection of a standard UV divergence of the free 
energy in the interacting boundary gauge theory. It is regulated by cutting off the 
radial coordinate integration at y — y c 1. It is straightforward to verify that the 
regularized Euclidean gravitational Lagrangian, Cf eg , is a total derivative, 


C E 

*^reg 


47 tG 5 
voi(n 3 
An Gh 


vol(H 3 


ryc d 

dy 

l dy 


4(1 ~y ) e 

y 2 


2 „-d 


4(1 -p) 2 e 

y 2 

a + 2 yad' 


2 „-d 


(a + 2 yad' 


ya 


yc 



(4.25) 


where in the second equality, using (4.23), we observe that the only contribution comes 
from the upper limit of integration. Regularized Lagrangian (4.25) has to be supplc- 
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merited with contributions coming from the familiar Gibbons-Hawking term, £q H , 


QH, _ 

bGH ~ 87 tG. 


5 JdMb 


= / dt E C, 


E 

GH i 


£ 


GH 


vol(h 2 3 
47 :Gk 


4(1 — y )e 

y 2 


—d 


( a{y - 4) - 2d'y(l - y)a + a'y{ 1 - y)) 


Vc 


(4.26) 


and the counterterm Lagrangian 10 , £ 


counter ? 


nE _ 

^counter 


dt e £ 


E 

counter ) 


£ s 

counter 


volQ 3 

471^5 




3 1 1 „ „ 3 a 2 

i + i Ri + 2 x +3Q '~2ta7: 


ln e«|-lx 2 ft-|x 1 )+lx 4 


3/c 


(4.27) 


where R 4 = R^He) is the Ricci scalar constructed from h E , and e c parameterizes 
conformal anomaly terms in terms of the gt B t E metric component, 




3 y 

2 ( 1 -!/)’ 



(4.28) 


The renormalized Lagrangian £® nom , finite in the limit t/ c —>- 0, is identified with the 
free energy T of the boundary gauge theory, 


T = £ 


E 

renom 


= lim ( £(L + Cr-!H + £? 

?/c—>0 


volh23 3 , 2 /"4 16 


- 1 + cT n -In 2 

4vrG 5 2 V 1>0 V 9 9 


Cl.o 


-3-3 in 2 ) +64p^ln2 


+ ^64piipio + g c 2,o + 32 p 10 — 4a 2 y 
lEli 1 + ~ ^(1 + 21n2 ){mtf + /32pf 0 + y(mf) 2 pi,o + jjc 2)0 - 4a 2)0 


(4.29) 


where in the second line we used the asymptotic expansion (4.9) and expressed the last 
line in terms of gauge theory variables using (3.4) and (4.13) and restoring the size t 
of the S 3 . Several comments are in order: 

■ For static gravitational solutions without Schwarzschild horizon (as discussed here), 
10 We keep only the counterterms relevant for the Rx S 3 background geometry of the gauge theory. 
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the free energy T must coincide with the energy E of the boundary stress-energy tensor. 
We explicitly verified that, indeed, 

T = E = E vacuum (m£) . (4.30) 

The latter is identified with the vacuum energy of Af = 2* gauge theory on S 3 . 

■ In a limit when all the (non-)normalizable coefficients vanish we recover the vacuum 
energy of the Af = 4 SYM (4.1). 

■ It is easy to extend discussion for general SO (4)-invariant non-equilibrium states of 
Af = 2* gauge theory — the final answer is as (4.29), except with {pi j0 , c 2 ,o , a 2 , 0 } now 
being functions of time. Note that 

^ °c ( 4|l6pi,ipi i0 (i) + ^c 2 , 0 (t) + 8p 2 lfi (t) -a 2 , 0 (t)J ^ = 0, (4.31) 

according to (4.12). That is, the boundary gauge theory energy conservation is enforced 
by the bulk momentum constraint (4.6). 


4.3 Vacuum states for rnf <C 1 


In preparation to the full numerical solution of (4.22), we discuss here its perturbative 
solution for p lvL <C 1. We introduce 


c = cosh(2A xi(y) + 0 ( A 3 )), p = e ^^iy)+o(^) ? 

a = 1 + A 2 a 2 (y) + 0 ( A 4 ), d = A 2 d 2 (y) + 0 { A 2 ), 


where A is a small parameter. Substituting (4.32) into (4.22) we find 

1+2/ / 3 


0 = Xi - 


0 = 0-2 


2/(1 -y) 
1 + 2 / 


xl + 


Xi 


% 2 (i -y) 

, 1 

a 2 + ~TT, - _ T Q '2 , 


2/(1 -y) ~ 2 / 2 (l- 2 /) 


0 = 4 - - IwC 1 - 9)(x',) 2 + -(xi ) 2 . 

2/(1 ~y) 3 y 


!\2 


0 = 4- ^2/(1 — 2/)(Xi) 


(4.32) 


(4.33) 


Solutions to (4.33) must satisfy boundary conditions corresponding to (4.9) and (4.23). 
We can solve equation for a 2 analytically, 


y In y 

«2 = Pi,1,(2) z -> 

1 — 2 / 


(4.34) 
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where Pi,i.( 2 ) is the non-normalizable integration coefficient. The remaining equations 
in (4.33) are solved with “shooting method” developed in [22], In particular, given the 
asymptotic expansions in the UV, y —> 0 + , 

xi =y 1/2 (i + y ^xi,o,(i) + \ ln y'j + °(y 2 ln y)^j , 

a 2 =^y + y 2 (a 2 , 0,(2) + ^ In y^j + 0 (y 3 In 2 y ), 
d-2=^y + y 2 +2xi,o,(l)-h^ln y'j + 0(y 3 ln 2 y), 


and in the IR, z —» 0+, 


Xi — Xo,(i) — > 


5 


«2 = (X 0 ,( 1 )) ^ - g* + 0 (z )J , 
d 2 =4,(2) ~~ ^(Xo^ 1 )) 2 ^ 2 + ^(4) j 


we hnd numerically, 


Ah, o,(i) 

1*2,0,(2) 

Xo,(i) 

4,(2) 

0.0568528 

-0.363452 

0.785398 

0.199266 


(4.35) 


(4.36) 


(4.37) 


To compare with the full numerical solution, we identify, to order 0( A 2 ), 

Pi,i = Pi,i,( 2 )A 2 , ci )0 = 2A 2 , pi,o = 0 , c 2 ,o = 4xi )0 ,(i)A 2 , 

a 2,0 = ®2,0,(2)A 2 , Po = 1 — Pl,l,(2)A 2 , Co = 1 + 2(Xo,(l)) 2 ^ , dp = 4(2)A 2 . 

(4.38) 

Note that AT = 2 supersymmetry in the UV at (9(A 2 ) leads to (see (4.14)) 


Pi.1,(2) - g 


From (4.29), 


e = 


Evacuum 

e m=a 

vacuum 


—1 + ^yrXi,o,(i) — 4a 2> o,(2) — -(1 + 2In 2)^ A" + 0 ( A 4 ) 

=1 + ^-Xi,o,(i) — a 2 ,o,( 2 ) — -(1 + 2 In 2)^ ( m £) 2 + 0((m£) 4 ). 


(4.39) 


(4.40) 
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C 2,0 Pl,0 







P 1,1 


Pl,l 


P 1,1 


Figure 1: Normalizable coefficients (4.24) as functions of pi i. The dashed lines repre¬ 
sent perturbative predictions (4.38) with (4.37). 

4.4 Gravitational solution and E vacuum for general ml 

Using the shooting method of [22], we solve (4.22) and determine the normalizable 
coefficients (4.24) as a function of mi = (12pi,i) 1//2 . The results of the computations 
for small values of pi_i are collected for numerical test in figure 1. The solid curves 
are obtained from numerical solution of full nonlinear equations (4.22), and the dashed 
lines represent perturbative prediction (4.38) with (4.37). 

In full nonlinear numerical analysis we constructed vacua for 0 < ml <8.5. The 
vacuum energy of the J\T = 2* gauge theory on S 3 relative to Af = 4 SYM Casimir 
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Figure 2: Vacuum energy of the AT = 2* gauge theory on S 3 relative to AT = 4 SYM 
Casimir energy, see (4.41). The vertical red line marks vanishing of e, see (4.42). 


energy is given by 

ffvacuum(m(.) 


e = 


E!f=± 


vacuum 


=l + AAt - |(1 + 21n2)(r«£) 2 ] 

+ ( 32p^ 0 + —(m£) 2 pio + -C2,o — ^ a 2 ,i 
It is presented in figure 2. The vertical red line indicates the mass scale mo£, 


(4,41) 


e(m 0 l) = 0 =4 m 0 e ss 0.87031, (4.42) 

at which the vacuum energy of the Af = 2* gauge theory vanishes and becomes negative 
for even larger value of rri£. 


5 Stability of J\f = 2* vacuum states within BEFP 

In the previous section we constructed gravitational solutions within PW effective 
action, identified as vacua of the Af = 2* gauge theory on S 3 . While the complete 
stability analysis of these solutions is beyond the scope of this paper, here we would 
like to analyze their stability within BEFP effective action. 

Effective action describing the fluctuations of an arbitrary PW static solution within 
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BEFP has been constructed in [12], 

SC = Cbefp — C PW F F 8Cy , 

SC 2 = -(1 + c) 2 (9X 2 ) 2 - i±£ ((c 2 + c)pf - 4(1 + c)p‘ /3 + dM) (X 2 ) 2 , 

W = -(1 + c) 2 (S.\V) 2 - ((c 2 - 1 )pf - 4(1 + c)pj /3 + (AV) 2 , 

(5.1) 

where p 6 = p 6 and X v = (X 3 ,X 4 , X 5 ) (see section 3 for more details). Note that 5C 
is SU{2)y invariant; as a result it is enough to consider a spectrum of only one of Xy 
components. In what follows we choose the latter to be X 3 . 

Introducing 

X 2 = e-^F 2 (y)n s (S 3 ), X 3 (t,y) = e~^F 3 (y)n s (S 3 ), (5.2) 


where Q s (S 3 ) are S 3 Laplace-Beltrami operator eigenfunctions with eigenvalues s = 
1(1 + 2) for integer l, 


A s3 tt s (S 3 ) = -s tt s (S 3 ) = -1(1 + 2) O^S 3 ), 
we find from (5.1) the following equations of motion 


(5.3) 


0 =F" + F, 


, ( 2 cc' (c 2 — l)p 8 4cp 2 2y — 1 1 

F ~ . . F / \ F 


2 

f 2 


c + 1 


6 ay 


e 2d co 2 s 


Fy(l — y)a \ a l-y 


3 ay y(y - 1) a(y - 1) 3 ap 4 y 

+ F 2 


c') 2 + p 2 (p 6 c- 4) 


(5.4) 


0 =F” + F: 


* [ 2cc ' + (c 2 ~ l)d 8 _ fcp^ + 2y - 1 + 1 


(l-c 2 )(c + l) 4ap 2 (p — 1)/ ’ 

2 


3 

^3 


C + 1 


6 ay 


3 ay y(y - 1) a(y - 1) 3 ap A y 


e u 


Fy(l ~y)a\ a l-y 


F F 3 


+ 


p 2 (p 6 (c-l)-4) 


(1 — c 2 )(c+l) 4ap 2 (p — 1) 


(5.5) 


The radial wavefunctions F 2)3 must be regular at the origin, he., z —> 0 + , 

F 2 = z 1 ' 3 p (1 + 0(z)) , Ft = z 1 ' 2 p (1 + 0(z)) , (5.6) 

and normalizable as y —> 0 + , 

9 — lu 2 ' 


f- 2 =y 3/2 (i Fy(^~ \ c ho + + °(y 2 ln y)^j > 


F 3 = y ( 1 + y ( - + 4 


s 4 _ w 2 

+ 4pi ;1 - 2p lj0 - -Ci,o - 2pi i In y ) + 0(y 2 In y) 


(5.7) 
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Figure 3: Low energy states in the spectrum of BEFP fluctuations about PW vacua: 
{n, 1} = {(0,0); (0,1); (1,0)} (blue, red, green). See section 5. 


Note that we set the normalizable coefficient of F 2.3 in the UV to one. 

When both scalars of the PW flow are set to zero, (5.4)-(5.7) corresponds to fluc¬ 
tuations of gravitational modes dual to dimension-3 (for F 2 ) and dimension-2 (for F 3 ) 
operators of the AT = 4 SYM on S 3 . In this case the equations can be solved analyti¬ 
cally. We fold, 


F 2 % M l} =y 3/2 ( 1 - y) l/2 2 F 1 ^-n,3 + n + l-l + 2;l- y y 

w 2 S m =3 + 2 n + l, 


(5.8) 


Fi =2/(1 - y) l/2 2^1 f-n ,2 + n + l ;l + 2 ;l-y) , 

V / (5.9) 

n!i} =2 + 2n + /, 

where {n, 1} are non-negative integers. For supersymmetric PW flows (4.14) we have 
to resort to numerics. The results of the numerical analysis are presented in figure 3. 
We look at the states with {n,/} = {(0,0); (0,1); (1,0)} for both F 2 and F 3 radial 
functions. Over the range of parameters discussed, the embedding of PW flows within 
BEFP effective action is stable. 


6 Black hole spectrum in PW effective action 

We begin with the metric ansatz and the boundary conditions representing regular 
Schwarzschild black hole solutions in PW effective action with the S 3 horizon. We 
explain how the normalizable coefficients of the gravitational solution encode the ther¬ 
modynamic properties of the black holes: the temperature Tbh, the energy Ebh, the 
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entropy S BB and the free energy E BB - We define the size i BH of a black hole as 

3 


BH 


A 


horizon 


L 3 


( 6 . 1 ) 


We compute excitation energy A (i BH /L, (ml )), 


A (i BH /L,(mi)) 


E BH (i BH /L, mi) - E vacuum {mi) 
E{?=* 


( 6 . 2 ) 


as a function of i BB /L, but for select values of mi: 

■ perturbatively in mi, to order 0((mi) 2 )] 

■ for pi i = A(mi) 2 = {1, 1.5, 2, • • • 5, 5.5, 5.8} (the last value corresponds to the largest 
value of mi for which we computed E vacuum ); 

and present a strong numerical evidence that 


lim A(i BH /L , (mi)) = 0 . (6.3) 

(■bh/L^ 0 

Thus, we conclude that there is no gap in the spectrum of black holes in PW geometry; 
correspondingly, there is no gap in SO (4)-invariant equilibrium states of the J\f = 2* 
gauge theory on S 3 in the planar limit and for large ’t Hooft coupling, as there is no 
energy gap for generic S'0(4)-invariant excitations in this theory. 


6.1 Metric ansatz and the boundary conditions for black holes in PW 


Recall that the vacuum solutions of section 4 were obtained within metric ansatz (4.2), 

4 


dsl 


vacuum 


COS 2 X 


-ae 2d (dt) 2 + — X - + sin 2 x(dEl 3 ) 2 


= - —ae 


-2d 


(< dt ) 2 + 


( dy) 


+ (i — y){dEt 3) 


(6.4) 


%(1 - y)a 

where in the second line we recalled the radial coordinate y = cos 2 a:, y G [0,1]. 
Regularity at the origin (y 1_) required that the metric functions a and d remain 
finite and non-zero. Notice that the three-sphere shrinks to zero size in this limit. 

In close analogy to (6.4), to describe regular horizon black holes, we reparameterize 
the radial coordinate y —>■ ypy, with a constant 0 < yh < 1, while keeping y G [0,1]. 
We further require that a has a simple zero and d remains finite as y —> 1_: 


ds 2 5 


4 

B h yhy 

0 < y h < 1, 


-ae 


-2 d 


(dt)' 2 + 


Vh{dyf 


y e [0,1], 


%( 1 - yyh)a 
lim a = 0 , 

y-U_ 


+ (1 - yyh)(dEL 3 ) 


(6.5) 


lim a' = finite A 0 , 
2 /^ 1 - 


lim d = hnite . 
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Given (6.5), 


-4 horizon 167T 


2(1-a -) 372 


1 4 1 '' 3 

BH _ horizon 


3/2 

2/fc 


(27t 2 ) 1/3 — . (6.6) 


1/2 

3/a 


The equations of motion describing black holes (6.5) can be obtained from (4.22) with 
the simple change of variables 11 y —» yyh, 


0 = d 


c ( c 'f . /Ac - l)p 8 4cp 2 a(2yy h - 1) + yy h 


+ d 


c 2 — 1 
p 2 (c 2 — l)(p 6 c — 4) 
4(1 - yyh)y 2 a 


6 ay 


+ 


3ay ya(yy h - 1) 3 yap 4 


0 = P n 


(P') 2 , / f (c 2 - l)p 8 4cp 2 a(2yyt - 1) + yy h 


+ P 


+ 


6 oy 3ay ya(yy ft -1) 3yap 4 
1 — p 6 c 


0 = d 


P 

(c 2 — l)p 9 
12(1 - yyh)y 2 a 6 p 3 y 2 a(l - yy h ) ’ 

„ 2y(l-yy /l )(c') 2 8(1 - yy h )y(p') 2 


3(c 2 - 1) 


P^ 


/ / 2 4 A(p') 2 2(c') 2 

0 — a — {y — yy h )a (— + J + ,(1 _ m) 


(Wh ~ 2 )a + yyh (c 2 - l)p 8 - 8p 2 c 


6y 


3yp 4 


(6.7) 
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We used the last two equations to algebraically eliminate a' and d! from the first two. 
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The boundary conditions in the UV, i.e., y —> 0 + , specify the asymptotic expansion 
a=l + y ~Ci,o + V 2 (a 2 ,o + Qci >0 (ci )0 + y h ) + 8 p 2 * + 16pi ; ipi )0 ^ hi y 


+ 8p\ x In y ) +0(y 3 \n 3 y ), 


d =y gCi j0 + y 2 


1. 1 - - 1. , 

2 C 2,0 — 3g C l,0 + 4 Pl,0 ~ g C l,0 yh + 2p 1;1 + 4pi )0 pi )1 


+ + ^ci )0 + 4p? ;1 + 8pi, 0 pi,i) hi y + 4 di,i ln2 y) + °(y 3 ln3 y) > 


p =1 + p (pi,o + Pi,i ln y) + y ( —c 10 + p lfi y h - 3pi,iCi )0 + 6p x x 


( 6 . 8 ) 


— 4pi, 0 pi,i + -ci,opi,o + -Pi,o + ^gPipCpo + Pi,iyh ~ 4p 1;1 

+ 3pi, 0 pi,i^ In y + ^p 2 ^ In 2 + £>(p 3 In 3 y ), 

c =1 + y ci, 0 + y 2 ^c 2 , 0 + + jjc 2 Inp^J + 0(p 3 In 2 y). 

In (6.8) the non-normalizable coefficients pi i and c 10 are related to corresponding 
coefficients of the vacuum solution as 


Pip — yhPip , ci,o — yhCifl , 


(6.9) 


to be further matched with the mass parameters {mbprrif} of the dual gauge theory 
as in (4.15). The rest of the coefficients in (6.8) are normalizable. The asymptotic 
expansion in the IR, i.e., as z = (1 — y) —> 0 + is different from the one in (4.23) - 
here it reflects the presence of a regular horizon (see (6.5)), 




- (eg) 2 ) U) s + 8 cg(pg) 2 + 


4 

W + 


6 yh \ 

l-Vh) 


d = d h 0 + 0(z ), 

P = Po + ®( z ) , 
c = c h 0 + 0(z ). 


+ 0(z 2 ), 


( 6 . 10 ) 


The full set of the non-normalizable coefficients is 


{«2,o , Pi,o , c 2 ,o, Po , Co , dg} . (6-11) 

Note that we have the correct number of non-normalizable coefficients to uniquely 
specify a solution of two second-order and two first-order ODEs given a choice of (6.9). 


24 



6.1.1 Perturbative black holes solutions 


As in section 4.3, we can construct solutions to (6.7)-(6.10) perturbatively in m£ to 
order 0((m£) 2 ). 

We introduce 


c = cosh( 2 Axi(z/) + 0 (A 3 )), p = e ^ 2 “2(?/)+o(A 4 ) ? 

a — (1 y)i l + V(l - Vh ^ + A 2 d 2 {y) + 0( A 4 ), d = \ 2 d 2 (y) + 0(\ 2 ), 
i - mjh 

where A is a small parameter. Substituting (4.32) into (4.22) we find 

o = v" - *1 (\ + + vi 1 ~ yh)(( 2 - y)yyh - z) \ _%_ 

Xl y(i - y) v (1 - Wh){ 1 +1/(1 - y fc )) 7 V(i - y )(i + y (i - yfc )) ’ 

n _ _ a 2 A , , y(i -yh)((2-y)m-2) \ _«2_ 

2 2/(1 -2/) V (1 - m)( 1 +1/(1 - y fc )) y 1/ 2 (1 - y)(i + y(i ^ y*)) ’ 

0 = al 2 - y — V -^ L -d 2 - |y(l - y)(l + 1/(1 - yh))(x’i ) 2 + ~(xi ) 2 , 

1/(1 - Wh) 3 y 

0 = 4 - |y(i -m)(xi) 2 - 

(6.13) 


For the asymptotic expansions we have: 
■as y-> 0+, 


Xi =y 1/2 (l + y (xi,o,(i) + In y ) + 0 ( 2/ 2 In 1/)) , 

«2 =Pi,i,(2) ^(tti,o,(2) + Iny) y + C>(i / 2 In y)^ , 

02 =^y + y 2 ^02,0,(2) + In y^ + 0 (y 3 In 2 y ), 

d 2 =ly + y 2 + 2 xi,o,(i) + y lny) + 0 (y 3 In 2 y ), 


(6.14) 


■ as z —> 0+ 

= ^ (1) (* ' + ° {z2) ) ' 

«2 =Pl,l,(2) ^«0,(2) ^1 “ ^ + C> (^ 2 )^ ’ 

02 =2(Xo,(i)) 22: + C , (^ 2 ), 

d 2 =do,(2) - 2(2 - + ^(^ 2 ) • 
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Equations (6.13)-(6.14) have to be solved numerically for different values of yh- 
To compare with the full numerical solution, we identify, to order 0( A 2 ), 

P 1,1 = Pi,i,( 2)A 2 , Ci,0 = 2A" , pi,o = Pi,i.(2)di,o,(2)A 2 , 62,0 = 4xi ; o,(i)A 2 , 

^2,0 — Vh — 1 + h2,0,(2)^ 2 , Po — 1 + Pl,l,(2)®0,(2)A 2 :) 

Co = 1 + 2(xo i(1) ) 2 A 2 , d h 0 = d h 0X2) A 2 . 

(6.16) 


Note that Af = 2 supersymmetry in the UV at 0( A 2 ) leads to (see (4.14)) 

1 

Pi,1,(2) - - • 


(6.17) 


6.2 Thermodynamic properties of black holes in PW 


Requiring that there is no conical singularity in the analytical continuation t —> He of 
the metric (6.5) as y —» 1_ we compute the Hawking temperature Tbh of the black 
hole using (6.10), 


T, 


_jh 

e 0 


BH 


- y h y / 2 
, 4(1 — yh) 


(i - y h )( 1 - WYMT + KY - y h mr + 6 y h 


C P h o ) 4 


(6.18) 


The Bekenstein-Hawking entropy of the black hole is given by 

-*4 horizon 4:71 (1 P/i) ^ 


S’ru — 


4G, 


Gr 


3/2 

2/h 


(6.19) 


The free energy Tbh can be computed following holographic renormalization procedure 
discussed in section 4.2. We find 


Tbh = 


3n 
4G7 


1 + 


1,0 


Vt V9 9 


16 In 2 + ^ In y h ) + Cl ’° 


9 




4 8, 4, 

.--- ln 2 +-lnp /l 


+ 32—y- (2 In 2 — In p/J 4—64pi ipi o + -C 2 o + 32p 2 0 — 4a2,o 
Vh Vh l 3 

(1 - y h )ne~ d o 

3 y 2 h G 5 


(1 - y h ){ 1 - (7o) 2 )(Po) 8 + 8c&(l - p,)(p^) 2 + 6 y h + 


/i\2 


4(1 — yh) 

(p h o) 4 


( 6 . 20 ) 
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The contribution in the last line in (6.20) comes from the lower limit of integration 
of the bulk contribution to the regularized free energy, (4.25); it equals precisely to 
(— SbhTbh )• Computing the holographic stress-energy tensor, as described in [10] we 
find 


Erh 


i,o (4 16 


V9 9 


ego / 4 8 


1 + -± n - IT ln2 + n ln Vh + — -o - o ln2 + - In y h 


3 3 


+ 32—j- (21n2 — In y h ) + — < 64p 11 p 10 + -£ 2,0 + 32/5^ 0 — 45.2,0 r ) 

Vh Vh.K ^ J / Ig 21) 

^ 1+ (fl.l (1+2ln2 . lnyh)(mef 

+ ^2“ ^32pj o + y(m^)Vpi,o + -c 2 ,o — 4a2,o, 

where in the last line we expressed the energy in terms of the dual gauge theory 
variables using (6.9) and (4.13). Notice that the basic thermodynamic relation, 


Ebh — Ebh ~ SbhTbh , 


( 6 . 22 ) 


is satisfied automatically. 

Using (6.16), from (6.21) we have 

E bh 1 , 4(1-j/ h ) / 8 . 


Ebh 

L 


=1 + 


( 812 \ 

X— Xi,o,(i)-02,0,(2) - -(1 + 2 In 2 - lnj/ fc ) (mfh 

3y/i Vh o J 


C>((m£) 4 ). 


(6.23) 


6.3 A (£ BH /L,{m£)) 

We are now ready to present results for A {£ B h/ L , (m£)) as defined by (6.2). 
To order 0((m£) 2 ), using (4.40) and (6.23), we find 

a 4 (1 -Vh) , a 2 i n\\\ 


+ A 2 (mi) 2 + 0({m£y ), 


A ! =A,( Sl )d(hffl- Iwl) 
3 \ Vh 


a 2,0,(2) \ 2 

- a 2,0,(2) 1 + X “W/i ' 

Vh J 3 


(6.24) 


Results of numerical computations of A 2 are presented in figure 4. A solid line repre¬ 
sents the data points, and the red dotted line is the best quadratic fit using the first 
10% of data points: 


A 2 =-0.0269118 


(6.25) 
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V-bhIL 



Figure 4: Solid line represents A 2 as defined in (6.24). The dotted red line represents 
the best quadratic fit to the first 10% of data points, see (6.25). 


Our numerical results present a strong evidence that 

lim A 2 = 0 , (6.26) 

(■bh/L^ 0 

as a result, we see that A vanishes in this limit to order 0((m£) 2 ). 

Using (4.29) and (6.21) we compute A for p 1;1 = ^(mi) 2 = {1,1.5, 2, • • - 5, 5.5, 5.8}. 
The results are presented in the left panel of figure 5 (the top-to-bottom blue curves 
correspond to p \^ variation 1 —> 5.8). The green curve represents A (mt = 0): 


2 4/3 / t RH 

A(m£ = 0 ) = - J73 (-^ 


+ 


2 2 / 3 ft 


BH 


(6.27) 


L J ' 7T 8 / 3 V L 

The right panel represents A for the largest value of mt computed: mt = 8.34266, with 
the red dotted line indicating the best quadratic fit to the first 10% of data points: 


A (mt = 8.34266) 


fit 


= 0.339765 f 


(6.28) 


Note that for mt = 8.34266, e = —243.785, implying that for the smallest size black 
hole studied, tsH/L = 0.0855056, 

Ebh — E v 


J vacuum 


= 1.04285 x 10 


-5 


(6.29) 


We conclude that numerical results strongly suggest (6.3). 
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A {mi = 8.34266) 
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Figure 5: Left panel: Black hole mass gap relative to E^£ um , see (6.2), as a function 
°f £bh/ L for select values of ml. The green curve represents A {ml = 0). Right panel: 
A for the largest value of ml computed, mi = 8.34266; the dotted red line represents 
the best quadratic fit to the first 10% of data points, see (6.28). 
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